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Abstract 

Simulations of nematic-isotropic transition of liquid crystals in two dimensions 
are performed using an 0{2) vector model characterised by non linear nearest neigh- 
bour spin interaction governed by the fourth Legendre polynomial P4. The system 
is studied through standard Finite-Size Scaling and conformal rescaling of density 
profiles of correlation functions. A topological transition between a paramagnetic 
phase at high temperature and a critical phase at low temperature is observed. The 
low temperature limit is discussed in the spin wave approximation and confirms the 
numerical results. 

PACS: 05.50.+q Lattice theory and statistics (Ising, Potts, etc.); 05.70.Jk Crit- 
ical point phenomena; 64.60.Fr Equilibrium properties near critical points, critical 
exponents; 64.70.Md Transitions in liquid crystals. 

1 Introduction 

The molecules of liquid crystals may often be described by long, neutral rigid rods which 
interact through electrostatic dipolar or higher order multi-polar interactions. This is at 
the origin of the natural introduction of Legendre polynomials for the description of the 
orientational transition between a disordered, isotropic high temperature phase and an 



ordered nematic phase at lower temperature [1]. When such a material is cooled at even 
lower temperatures, other ordered phases may be encountered, the description of which 
requires more realistic potentials (see e.g. Ref. [2]). 

Lattice models of nematic-isotropic transitions capture the essentials of the above 
description. The molecules are represented by n— component unit vectors <x r , hereafter 
called "spins", located on the sites r of a simple hyper-cubic lattice. The interaction 
between molecules is restricted to the nearest neighbours (r, r'), so that the radial depen- 
dence is kept constant and the angular dependence enters through a k— th order Legendre 
polynomial 1 , Pk(cosa r y), where a r y = (cr r ,cr r /) is the angle between vectors cr r and 
cr r /. The intensity of the interaction energy is measured through a parameter e. The 
Hamiltonian of a lattice liquid crystal is thus given by 



and the relevant parameters for the investigation of the phase transition are the space 
dimension D, the "spin" dimensionality or equivalently, its symmetry O(n), and the 
"symmetry" k of the interaction 

Usual XY and Heisenberg models correspond, within this description, to k — 1 and 
n = 2 and 3, respectively. In two dimensions (D = 2), these models exhibit quite different 
behaviours. In the case of the XY model (0(2) abelian symmetry), a topological tran- 
sition described by Berezinskh, and Kosterlitz and Thouless [3-5] takes place, governed 
by the condensation of vortices. This is not prevented by the Mermin- Wagner- Hohenberg 
theorem [6-8] which states that spin models with continuous symmetry cannot have any 
long range ordered phase. To leading order in the high-temperature expansion, one gets 
for the correlation function an exponential decay (<r ri -cr r2 ) ~ K^ ri ^ r ' 2 ^ (/f — e/A^T), while 
at low temperatures in the harmonic approximation, the Hamiltonian becomes quadratic 
and leads to the Gaussian model which implies [9] that (<r ri ■ cr r2 ) ~ |ri — y 2 \^ 1 ^ 2 ' kK , i-e. 
a temperature-dependent spin-spin critical exponent t]xy(T) = T — > 0. The low 

temperature (LT) phase of the XY model is a quasi-long-range ordered phase (QLRO) 
with vanishing magnetisation M 2 (T) = lim| ri _ r2 |^ 00 (cr ri ■ cr r2 ) = 0, or a critical phase. 
The Heisenberg model (0(3), no n- abelian symmetry) on the other hand has no transition 
at any finite temperature (asymptotic freedom) [10-13]. This difference with XY model 
is at first surprising, since the low temperature limit of the Heisenberg model is essen- 
tially described by a similar spin wave approximation (SWA): the longitudinal modes 
of the O(n) model are frozen and only the transverse modes are activated, leading es- 
sentially to two Gaussian models. This apparent contradiction between the asymptotic 
freedom of the Heisenberg model and the topological transition at finite temperature for 
the XY model finds its origins in the stability of topological defects in the latter case, 
while the 'third spin dimension' makes the vortices unstable at any temperature in the 
former model. Further, the question of the accessibility of the thermodynamic limit is 
worth studying. Berezinskh and Blank noticed long time ago that a really large, but finite 
XY system always possesses a non-zero magnetisation [14,15]. The size being limited, 

1 Even order Legendre polynomials guarantee the local Zi symmetry er r — > er r . 
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1 ri — r2 1 becomes at most as large as the linear size L and a finite order parameter follows, 
M L (T) ~ L-^ (r) . More precisely, M L (T) = 0(1) as long as L « e 2/r?(T) a condition 
which can be fulfilled for any L by considering small enough temperatures. With this 
result in mind, we expect for the model considered hereafter that a spin wave solution 
will be found at low enough temperature and thus we search for an algebraic decay of 
the spin-spin correlation function from an effectively non-zero order parameter, although 
there is at most QLRO at low temperatures. 

Changing the value of k in Eq. modifies the symmetry of the spin-spin interaction 
and gives rise to new features 2 . When k increases, one may indeed expect a qualitative 
change in the nature of the transition, like in the case of discrete spin symmetries (Potts 
model) [18, 19]. The value k = 2 was intensively studied. It still corresponds to the XY 
model for 0(2) spin symmetry, while it leads to the RP 2 or Lebwohl-Lasher model [20] for 
3— component spin vectors. The nature of the transition in this latter case is still under 
discussion [21-25], but a recent study reported new evidences, extremely convincing, in 
favour of a topological transition [26]. The transition is driven by topologically stable 
point defects known as | disclination points. Considering still larger values of k, there 
is a proof of asymptotic freedom in the large— n limit, for values of k (in the interaction 
term (1 + cos6') fc ) which do not exceed a critical k c ~ 4.537... [23]. This is again discussed 
in a recent preprint from Caracciolo et al. [27]. Above this value the transition becomes of 
first order, a result which does not violate Mermin-Wagner-Hohenberg theorem, since the 
correlation length is finite at the transition. For finite value of n, the question of the nature 
of the transition at high k is still a challenging problem. In the context of orientational 
transitions in liquid crystals, Legendre polynomials rather than cos fc 6 interactions are 
introduced, and we are led to the Hamiltonian of Eq. (0). The 3- vector model with P4 
interactions was already considered in Refs. [28,29] where convincing evidence for a first 
order transition was reported. 

In this paper, we study the behaviour of an abelian spin model, namely 0(2) rotation 
group with P 4 — like spin interactions. The model will be referred to as P4 0(2) for 
simplicity. We are mainly interested in the low temperature properties of the model 
where comparisons with analytic predictions are possible due to the simplicity of the 
SW approximation. The techniques used combine temperature analysis, FSS (Finite-Size 
Scaling) and conformal techniques (Finite-Shape Scaling - FShS - to plagiarize the famous 
acronym) . 

2 Definition of the model and of the observables 

In Refs. [28,29], the following Hamiltonian Hp 4 0(3) = — e^( rr ') P^i^r ' °v) was consid- 
ered, where cr r = (a* , cr^ , a*) , \er r \ = 1 and P±(x) = |(35x 4 — 30x 2 + 3). For 2— component 
vectors, in the completely disordered phase (cos 2 9) = | and (cos 4 ^) = |. In order to 
keep the same symmetry in the interaction than in the P4 0(3) model, but to normalise 
it between and 1 in the limits of completely disordered and completely ordered phases 

2 In the same spirit, the case of symmetry-breaking magnetic fields hk cosk8 added to the XY model 
and changing the phase diagram was investigated by Jose et al [16, 17]. 
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respectively, we modify slightly the Hamiltonian to include pair interactions of the type 
Qa(x) = AP±{x) + const = ^(35x 4 — 30x 2 + ^). The corresponding Hamiltonian is thus 
defined by 

Hp A 0(2) = ~ e ^2 < 5 4 ( <Tr ' °"r')> ( 2 ) 

(r,r'> 

with now cr r = (cr*, <j^), |cr r | = 1. A qualitative description of the transition is provided by 
the temperature behaviour of the energy density, the specific heat, the order parameters 
and the corresponding susceptibilities. The internal energy is defined from the thermal 
average of the Hamiltonian density, 

u P4 o(2)(T) = (DL D )- 1 (H P4 o(2 ) ) (3) 
and the specific heat follows from fluctuation dissipation theorem, 

L D T 2 C V (T) = ((H P4 o (2 )) 2 ) - (H P4 o (2 )) 2 . (4) 

Brackets denote the thermal average. The definition of the scalar order parameter (some- 
times called nematisation) is deduced from the local second-rank order parameter tensor, 

q^( r ) = - h af >. (5) 

After space average, the traceless tensor L~ D J2 r <? Q/3 ( r ) admits two opposite eigenvalues 
±^rj corresponding to eigenvectors n + and n_. The order parameter density is defined 
after thermal averaging by 

Q2(T) = (v). (6) 

This quantity has the same physical content but is more stable numerically than a direct 
estimation of (2(<r r ■ n + ) 2 — 1). Another order parameter may be defined simply by 
inspection of the structure of the Hamiltonian, 
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(T)=L- B (^g 4 ( CTr .n + )). (7) 



The associated susceptibilities are defined by the fluctuations of the order parameter 
densities, e.g. 

4L D 

X q2 {T) = Y j{{ri 2 )-{ri) 2 )- (8) 



3 Thermal behaviour 

In this section, we illustrate the behaviour of the various thermodynamic quantities as 
the temperature varies. It gives a first idea of the nature of the transition. 

The simulations are performed using a standard Wolff algorithm suited to the expres- 
sion of the nearest neighbour interaction [29,30]. The spins are located on the vertices of 
a simple square lattice of size L 2 with periodic boundary conditions in the two directions. 
We use L = 24, 32, 48, 64 and 128 with 10 6 equilibrium steps (measured as the number 
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of flipped Wolff clusters) and 10 6 Monte Carlo steps (MCS) for the evaluation of thermal 
averages. The autocorrelation time (at k B T/e = 0.2, L = 16) is of order of 30 MCS, 
hence the numbers of iterations that we used correspond roughly to 3.10 4 independent 
measurements for the smallest size and is still safe at L = 128. 
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Figure 1: Energy, specific heat, order parameter q 2 and corresponding susceptibility x<j 2 > 
order parameter g 4 and corresponding susceptibility Xqa vs T for the P4 0(2) model. The 
full lines are only guides for the eyes. The values of k B and e have been fixed to unity. 



We deliberately did not use histogram reweighting which would require a huge amount 
of simulation time to be reliable. We have adopted a different strategy here, producing 
more simulations with less iterations, but with a better control of the errors than at the 
ends of the histograms. The temperature dependence of thermodynamic quantities is 
plotted in Fig. ^ for different system sizes. The behaviour of the energy density clearly 
displays a difference between the regimes of low and high temperatures. This is the 
signature of a transition and this naive statement is corroborated by the behaviours 
of the other physical quantities, the specific heat C v , the order parameters q 2 and q^, 
and the corresponding susceptibilities. The specific heat close to the maximum does 
not seem to increase substantially with the system size. This might be the sign of an 
essential singularity 3 around a temperature ksTc/e ~ 0.70 — 0.75. From the behaviour 

3 Or at least the sign of a non diverging specific heat with non-positive exponent a. 
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of the order parameters, we may suspect a smooth transition, since there is no sharp 
jump. The susceptibilities display a non conventional behaviour at low temperature, 
increasing with the system size, which indicates a likely topological transition with a 
critical low temperature phase where the susceptibility diverges at any temperature (note 
the logarithmic scale for vertical axis). 

The probability distributions of the energy density and the order parameter is also 
instructive. Both quantities are shown in Fig. |21 The distributions have a simple shape 
with single peaks, and this is still true for temperatures below the transition, a result 
which suggests a continuous transition. We note that due to the finite-size of the system, 
the order parameter is finite. 




Figure 2: Energy density and order parameter probability distributions (arbitrarily nor- 
malized) for the P4 0(2) model at a temperature close to the transition . The system size 
is L = 48 and the temperatures are given in units suche that feg/e = 1. The solid lines 
correspond to the effective transition. 



4 Finite-Size Scaling 

In this section we investigate the properties of the low temperature phase using Finite- 
Size Scaling technique. In the critical low temperature phase of a model which displays 
a topological transition (the paradigmatic XY model serves as a guide), the physical 
quantities behave like at criticality for a second-order phase transition, with power law 
behaviours of the system size. The difference is that in the critical phase, the critical 
exponents depend on the temperature and for any temperature below the transition one 
has e.g. 

q 2 {T) ~ L-!^( T ) (9) 
Xq2 {T)-L 2 ^ T \ (10) 
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Here i] Q2 (T) denotes the correlation function critical exponent, defined by 

(Q 2 (cos(# ri - r2 ))> ~ In - r 2 \-^ T \ (11) 




Figure 3: FSS behaviour of the order parameter densities q 2 (T) and (T) and the corre- 
sponding susceptibilities for the P 4 0(2) model at temperatures ksT/e = 0.1, 0.2 and 0.4 
below the transition temperature. Temperature-dependent exponents are confirmed and 
the full lines show the expected SW result which nicely fits the data at low temperatures. 

In Fig. |3]the FSS behaviour of the order parameter densities and of the corresponding 
susceptibilities is shown on log-log scales vs the system size L (we used L = 16, 32, 48, 64, 
128 and 256) for different values of T below the expected transition. The linear behaviour 
on this scale indicates power laws, and the different slopes at different temperatures is 
the result of exponents which depend on the value of the temperature. Starting from a 
fully ordered system at T = 0, at low temperature we can expect a small disorientation 
of the molecules and the spin- wave approximation becomes correct. For 0(2) model with 
nearest neighbour interactions described by arbitrary polynomial in cosa rir /, one is led to 
an effective harmonic Hamiltonian |e r A l(9 r — 9 r /) 2 . It yields power-law correlations, 

2 

(cosm(# ri - r J) = e'^ {i9 ^- 9 ^ 2) ~ |n - r 2 p^ (12) 
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with a spin-wave decay exponent given by 



n s J = -vMT) = (13) 

This expression, if confirmed numerically, will support the presence of a quasi-long-range 
ordered, scale- invariant phase at low temperatures. The exponent t] q2 {T) (m = 2) is 
accessible through q 2 (T) (or the corresponding susceptibility) via equations Q and (fTUj). 
We note that it is also accessible through the the FSS behaviour of 94 (T), since the leading 
behaviour of {Q±{6 Tl —9 r2 )) is still governed by (cos 2(9 Tl —0 T2 )). In FigEl one may extract 
numerical values of 77 exponent. The values rj q2 (T) = 0.0056, 0.0117, 0.0182, and 0.0256 
are obtained at kpT/e = 0.1, 0.2, 0.3, and 0.4, respectively. Using the small angle limit 

Hp 4 o(2) = - e S(r,r'> <54(cOSa r>r ~ |e 5^<r,r') TT^r ~ ^r') 2 ' We naVe 1 = TT wnicn yields 

^fSa = 1 iff • At ^ T / e = ai > °- 2 > °- 3 > and °- 4 > we s et - °- 006 > °- 011 > °- 018 > 

and 0.022 which confirm the numerical values, since the SW exponent is always a lower 
bound for the exact exponent rj q2 (T). The comparison is made visible in the figure where 
the SW values are plotted in full lines while the symbols represent the numerical data. 
As expected, the lower the temperature, the better the SW approximation. 



5 Finite-Shape Scaling 

Inspired by the results that we obtained for the Lebwohl-Lasher model [24] or the XY 
model [31-33], we will now produce a complementary study using a rescaling of the density 
profiles. For that purpose, we assume the existence of a critical phase at low temperatures 
as suggested by the temperature dependence and FSS results. If the assumption is revealed 
incorrect, we will be led to some inconsistency. 

The existence of a scale-invariant low temperature critical phase leads to conformally 
covariant density profiles or correlation functions at any temperature below the transition 
T c . It is then advantageous to deduce the functional expression of the correlation func- 
tions or density profiles in a restricted geometry adapted to numerical simulations from 
a conformal mapping w(z): 

G{ Wl ,w 2 ) = \w'{z 1 )\-^\w'{z 2 )\-^G(z 1 ,z 2 ) (14) 

Here, w labels the lattice sites in the transformed geometry (the one where the computa- 
tions are really performed), z is the corresponding point in the original one (usually the 
infinite plane where the two-point correlations take the standard power-law expression 
G(zi,z 2 ) ~ \zi — z 2 \~ r,a ), and x a = ^r] a is the scaling dimension associated to the scaling 
field under consideration. The interest of such an approach lies in the full inclusion of 
the changes due to shape effects in the functional expression and we may copy the termi- 
nology FSS and call it Finite-Shape Scaling. Rather than two-point correlation functions, 
it is even more convenient to work with density profiles m(w) in a finite system with 
symmetry breaking fields along some surfaces in order to induce a non-vanishing local 
order parameter in the bulk. The density m(w) will be q 2 (r) = (Q 2 (o'r • hsA)) or higher 
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Figure 4: Finite-Shape Scaling behaviour (L = 64) of the order parameter density profile 
52 (w) for various temperatures kgT/e below the transition temperature and one temper- 
ature above. The transition takes place around the value 0.70. The insert shows (on a 
semi-log scale) the exponential decay of correlation function (below the transition) along 
a torus as explained in the text. 



rank nematisation 94 (r) = (Qi{cr r ■ haA.))- in the case of a square lattice A of size L x L, 
with fixed boundary conditions along the four edges <9A, one expects 
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where w stands for lattice site r. This expression easily follows from the expression of the 
order parameter profile decaying in the upper half-plane from a distant surface of spins 
constantly fixed in a given direction, m(z) ~ y~ x ° , and from the conformal transformation 
of the upper half-plane z = x + iy, (0 < y < 00), inside a square w = u + iv of size L x L, 
(-L/2 < u < L/2,0 < v < L), with open boundary conditions along the four edges, 
realized by a Schwarz-Christoffel transformation 



w(z) = —F(z,k), 



sn 



(16) 



Here, F(z, k) is the elliptic integral of the first kind, sn(2Kw/L) the Jacobian elliptic 
sine, K = K(k) =F(1, k) the complete elliptic integral of the first kind, and the modulus 
k = 0.171573 depends on the aspect ratio of A (here 1). 
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The procedure is now to fit numerical data of the order parameter profile against 
expression (|15|). The result for the density profile of qi (w) is shown on a log- log scale in 
Fig. EJ At low temperatures, the resulting straight lines confirm the existence of a critical 
phase (a linear behaviour on this scale results from an algebraic decay in the original 
semi-infinite geometry). Above the deconfining transition (assuming that it is indeed the 
driving mechanism of the transition), the decay becomes faster, indicating a paramagnetic 
phase. The transition is approximately located at a temperature kBT c /e ~ 0.70 — 0.75. 
The scenario is eventually completely consistent with a BKT transition. Furthermore the 
7] exponent again follows from the scaling of the density profile in Eq. (|15|) which provides 
an alternate determination of this quantity. 

Another famous conformal mapping which has been applied to many two-dimensional 
critical systems is the logarithmic transformation w(z) = ^\nz = ^ In p + i-^-. It maps 
the infinite plane onto an infinitely long cylinder of perimeter L, and due to the one- 
dimensional character of this latter geometry, the correlation functions along the axis of 
the cylinder (let say in terms of the variable u = ^- In p) decay exponentially at criticality, 
G(ui, U2) ~ exp[— («2— The interesting result which makes this technique powerful 

is that the correlation length amplitude on the strip is universal and only determined by 
the corresponding rj exponent, £ = This relation is known as the gap-exponent relation 
in the context of quantum chains in 1 + 1 dimensions. It was conjectured by several 
authors [34-36] before Cardy proved it [37]. Using MC simulations we cannot of course 
produce an infinitely long cylinder. It is however possible to perform simulations inside a 
rectangle L\ x L with L\~^> L and periodic boundary conditions in both space directions 
(we get a very long torus). Due to the exponential decay of the correlation length, it 
is not necessary to explore really long distances 112 — u\ (typically 112 — u\ < 10L). A 
finite torus of long perimeter 10 3 L thus only produces insignificant finite-size corrections 
to the gap-exponent relation. We performed the simulations at different temperatures 
in a system of size 10 x 10000 and extracted the correlation function exponent from the 
linear behaviour 

1X11 

ln(<3 2 [cos(0 U2 - 9 Ul \) = const ~j~{ u<1 ~ Ul )- ( 17 ) 

This linear behaviour (in terms of the variable u<i — U\) is shown in the insert of Fig. 0] 
where for simplicity (Q2[cos(9 U2 — 9 Ul ]) is denoted G 2 (u2 — Ui). It is impossible to apply 
this technique up to the transition temperature, since the strip system contains quite a 
large number of spins (10 5 while simulations in the square geometry are performed up to 
typically 10 4 spins) and the autocorrelation time increases too fast. 

6 Behaviour at the deconfining transition 

Not only the low temperature behaviour is interesting. The value of the r] exponent at 
the BKT transition where some deconfining mechanism should lead to the proliferation 
of unbinded topological defects is also of primary interest. For that purpose, an accurate 
value of the transition temperature is needed. We performed a study of the crossing of £7 4 
Binder cumulant for very large statistics (30 x 10 6 MCS) and large system sizes (squares 
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results shown in Fig. |3] 




Figure 5: Left: Crossing of the Binder parameter at the deconfining transition at a 
temperature k B T BKT /e = 0.7226. Right: Size dependence of the correlation function 
exponent at Tbkt {L = 16, 32, 64, 128) and extrapolation to the thermodynamic limit. 



Then this temperature is used to perform Finite-Shape Scaling using the technique 
already employed for the investigation of the low temperature phase, namely the alge- 
braic decay of density profiles inside a square with fixed bounday conditions. These new 
simulations are really time-consuming, since the autocorrelation time increases in the low 
temperature phase when T evolves towards the deconfining transition and a rather large 
number of Monte Carlo steps is needed to get a satisfying number of independent mea- 
surements. For sizes L = 16, 32, 64, we used 10 6 MCS for thermalization and 30. 10 6 
for measurements, while "only" 20. 10 6 for the largest size 128. For technical reasons, the 
exponential decay of two-point correlation functions along the torus cannot be applied at 
the BKT transition, since as already mentioned the system size being quite larger than 
in a square geometry, the number of MC iterations required is by far too large. In Fig. |5] 
we plot the "effective" exponent r) e g(L) measured at Tbkt for different system sizes as a 
function of the inverse size. An estimate of the thermodynamic limit value (L — > oo) can 
be made using a polynomial fit (the results of quadratic and cubic fits are respectively 
0.118 and 0.122). It is safer to keep the three largest sizes available, L = 32, 64, and 
128, for which a linear dependence of r] c ^{L) with L~ l is observed. Taking into account 
the error bars, crossing the extreme straight lines leads to the following value for the 
correlation function exponent at the deconfining transition 

% 2 (Tbkt) = 0.122 ±0.007. (18) 

This value is essentially half the Kosterlitz value for the XY model. 
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7 Conclusion 



In Fig. El we plot as a function of the temperature the exponent t] q2 (T) measured after 
conformal rescaling of the density profile and correlation functions at different system 
sizes and the FSS determination which follows from Fig. 01 Together with the exponent 
determined numerically, we report the result of the spin-wave approximation, shown in 
dotted line. The larger the size of the system, the better the agreement. Similar results 
(not shown here) are measured in the case of the higher-order nematisation, q±(T). 
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Figure 6: Temperature variation of the correlation function exponent r] q2 (T) deduced 
from conformal rescaling (open symbols) and FSS (filled symbols). Open triangles and 
diamonds, which correspond the largest systems where we appplied FShS seem quite 
reliable. The dashed lines shows the result of the spin-wave approximation i] q2 (T) = 
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According to these results, the main outcome of the present work is the following: 

- P 4 0(2) model displays a BKT-like transition with QLRO in the LT phase where SWA 
nicely fits the nematisation temperature-dependent exponent rj q2 (T) when T — » 0. 

The LT phase of the model is thus disordered. Nevertheless, in what we will call a physical 
limit, i.e. a finite but quite large size where the system contains a macroscopic number 
of spins, it is useful to consider that the system is partially ordered, as Fig. |21 seems to 
indicate. 

These results find a partial interpretation through a naive comparison with clock model 
in 2D. Increasing the order of the interaction polynomial indeed increases the number of 
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deep wells which stabilise the relative orientation of neighbouring spins. One is thus led 
to a system which is quite similar to a planar clock model with a finite number of states, 
unless the fact that here we keep a continuous spin symmetry which prevents from any 
"magnetic" long-range order at finite temperature. The clock model is known to be in the 
Potts universality class when q = 3, but at q > 4, it displays a QLRO phase before conven- 
tional ordering at lower temperatures [16]. Combining these results with the requirements 
of Mermin-Wagner-Hohenberg theorem in the case of continuous spin symmetry gives a 
natural framework for the comprehension of our results for 2— component spin systems. 
Whatever the nearest-neighbour interaction (in P 1: P 2 or P 4 ) their behaviour seems to be 
always described by a BKT transition. The similar observation that a two-component ne- 
matic model renormalises in two dimensions towards the XY model was already reported 
in Ref. [38]. The transition is likely driven by a mechanism of condensation of defects, like 
in the XY model, but due to the local Z 2 symmetry not only usual vortices carrying a 
charge ±1 are stable, but also disclination points carrying charges ±1/2 should be stable. 
The role of these defects might be studied in a similar way than in the recent work of 
Dutta and Roy [26], by the comparison of the the transition in the pure model and in a 
modified version where a chemical potential is artificially introduced in order to control 
the presence of defects. 
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